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ABSTRACT 



While the population of main sequence debris discs is well constrained, little is 
known about debris discs around evolved stars. This paper provides a theoretical 
framework considering the effects of stellar evolution on debris discs, particularly the 
production and loss of dust within them. Here we repeat a steady state model fit to disc 
evolution statistics for main sequence A stars, this time using realistic grain optical 
properties, then evolve that population to consider its detectability at later epochs. 
Our model predicts that debris discs around giant stars are harder to detect than 
on the main sequence because radiation pressure is more effective at removing small 
dust around higher luminosity stars. Just 12% of first ascent giants within lOOpc are 
predicted to have discs detectable with Herschel at IQOfini. However this is subject 
to the uncertain effect of sublimation on the disc, which we propose can thus be 
constrained with such observations. Our model also finds that the rapid decline in 
stellar luminosity results in only very young white dwarfs having luminous discs. As 
such systems are on average at larger distances they are hard to detect, but we predict 
that the stellar parameters most likely to yield a disc detection are a white dwarf at 
200pc with cooling age of O.lMyr, in line with observations of the Helix Nebula. Our 
model does not predict close-in (<0.01AU) dust, as observed for some white dwarfs, 
however we find that stellar wind drag leaves significant mass (^1O~^M0), in bodies 
up to '-^lOm in diameter, inside the disc at the end of the AGB phase which may 
replenish these discs. 

Key words: 



1 INTRODUCTION 

The first dusty disc around a main sequence star was ob- 
served in 1984 around Vega (|Aumann et al.l 1 19841 ). Since 
then our knowledge of such systems has improved signifi- 
cantly, and it is now known that 32% of A stars exhibit 
excess emission in the infrared , over and above the stel- 
lar photosphere (|Su et al.l 120061 ). This is thermal emission 
from dust particles orbiting the star in a debris disc. De- 
bris discs are coUisionally dominated in that the smallest 
bodies in the system are continuously replenished by colli- 
sions between larger objects and are subsequently removed 
by radiation pressure. The long term evolution of such sys- 
tems can be mode l led by steady state collisional models 
llKriyov et al.ll2008l: IWvatt et all [20071 : iKenvon fc B^omle^l 
I2OO4I : [Kennedy fc Wvattll2010l ). Disks are depleted due to 
collisional erosion and are expected to show a slow decline 
in brightness . A decrease i n brightness with age is indeed 
observed i.e (|Su et al.ll2006l ) and can be well fitted with the 
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models of lWvatt et al.l ((20071), allowing such models to char- 
acterize the population of main sequence A stars debris discs 
reasonably accurately. These models assume that velocities 
in the disc are high enough that collisions are destructive. 
This require s that the disc is st i rred, for exmple by self- 
stirr ing (e.g. Kenvon fc BromlevI (I2OOJ)) or planet stirring 
(e.g. iMustiU fc WvattI (|2009t ). 

Dust is also seen around some post-main sequence stars. 
In some cases this dust can be a result of the evolution of 
the star, for example material emitted in the stellar wind 
form spher ical shells of du s t that are observed around AGB 
stars (e.g. lOlofsson et al.l (|2010l )) or even stable discs ob- 
serve d around post-AGB stars , possibly linked to binarity 
(e.g. Ivan Winckel et al.l ( 20091 )) . Infr ared excess observed 
arou nd giant stars, e.g. Ijural (| 19991 ). and the helix neb- 
ula (|Su et al.l 120071 ). on the other hand, has been inter- 
preted as a disc similar to debris discs on the main se- 
que nce (although alt ernative interpretations do exist see 
e.g. iKim et al.l (200l|)). Hot dust is also obser ved in small 
radii ( <0. 01AU) discs around whit e dwarfs, e.g. lFarihi et al.l 
(|2009l ) or Ivon Hippel et al.l (|2007h , again inferred to origi- 
nate from a debris disc. However, in contrast to main se- 
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quence debris discs, these discs cannot be in steady state 
since material at such small radii has a short lifetime. Rather 
models suggest that these discs are formed when an aster- 
oid approac hes close to the star where it is tidally disrupted 
l|jurall2008l ). 

There are not yet enough observations of discs around 
post-main sequence stars to fully understand the population 
and it is not clear how the few discs that have been discov- 
ered around post-main sequence stars relate to the progeni- 
tor population of debris discs on the main sequence. In this 
work we take advantage of the fact that the main sequence 
debris disc population around A stars is well characterised 
and extend the steady state coUisional evolution models to 
consider the changes to this known population during the 
star's evolution. In particular we consider its detectability 
on the post-main sequence and therefore whether the ob- 
served post-main sequence discs derive from this population 
and what future observations would be best suited to detect 
them. 

Previous work has looked at specific aspects of the ef- 
fects of stellar ev olution on astero ids or comets, su ch as 
stella r wind drag (IDong et al.luOlOl ') and sublimation (jJural 
|2004 I2OO8I '). iDong et all (|201Cl ) model the evolution of a 
planetesimal belt due to stellar mass loss and suggest that 
the captue of km-sized planetesimals into mean motion res- 
onance could explain systems such as the helix nebula. In 
this work a theoretical framework is built that incorporates 
the effect of collisions, radiation forces, the stellar wind, sub- 
limation and realistic optical properties of dust, during the 
star's evolution from the main sequence to the white dwarf 
phase, focussing on the observable properties of the belt. 
The dynamical efi'ects of stellar evolution, in particular stel- 
lar mass loss, on planetary systems will be considered in 
future work. 

This paper begins by discussing the evolution of the star 
in ^ In ^we revisit the initial population of debris discs 
around main sequence A stars this time rep roducing the 
fit to Spitzer observations at 24 and 70/im of IWvatt et al.l 
((20071) using a model that incorporates the optical properties 
of realistic grains rather than a black body approximation . 
The steady state coUisional models of IWvatt et al.l (|2007l ) 
are then extended to include post-main sequence stellar evo- 
lution. In 52] the changes to the properties of an individual 
disc, as the star evolves, are described, whilst ij5] discusses 
the implications for observations of the population of debris 
discs around evolved stars, focusing in particular on giant 
stars and white dwarfs and !j6] summarises the models dis- 
cussed. 



2 STELLAR EVOLUTION 

Knowledge of the evolution of the star is required in order to 
model the post-main sequence evolution of the debris disc. 
The two crucial parameters that determine the disc's evolu- 
tion and detectability are the stellar lu minosity and mass. In 
this work stellar evolution models from lHurlev et al.l (|200Cl ) 
are used. These models use analytic formulae and are de- 
signed for population synthesis. They are accurate to within 
5% of more detailed models. 

Here we consider stellar evolution models for stars of 
mass 1.67 to 3.8 Mq, these correspond to stars of spectral 
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Figure 1. The luminosity (L), mass (M), mass loss rate and 
temperature (T) evolution of a 2.9Mq star, Z=0.02, in solar val- 
ues. MS: main sequence (0-413Myr), RGB: red giant (415.9-418.7 
Myr), HB: horizontal branch (418. 7-521. 4)Myr, AGB: asymptotic 
giant branch (521. 4-525. 9)Myr, WD: white dwarf (525.9 Myr on- 
wards). 



type A9-B8, according to the models of lKurucj (|l979l ). al- 
though the stellar temperature varies along the main se- 
quence and a star of a given mass will be classified differ- 
ently depending on the point in its evolution at which it is 
observed. 

The evolution of A stars can be split up into vari- 
ous phases (see Fig. [T] for the evolution of an 2.9Mq star, 
equivalent to spectral type AG). Along the main sequence 
(MS) the star is stable and powered by core hydrogen burn- 
ing. There is a slight increase in luminosity, as the size of 
the hydrogen burning core increases, but mass loss is neg- 
ligible. Main sequence luminosities, for the spectral type 
stars considered, range from 7.5Lq (1.67Af0) to 1801/0 
(3.8Mq) and main sequence lifetimes from 200Myr (3.8Mq) 
to 2000Myr (I.67M0). Changes in stellar properties along 
th e main sequence ar e negligible and therefore are ignored 
in IWvatt et al.l (|2007l ) and !j3l Once hydrogen in the core is 
exhausted, hydrogen shell burning commences and the star 
swells to become a red giant (RGB). Its luminosity increases 
by several orders of magnitude. The RGB phase lasts for 
~10Myr, until core temperatures are hot enough that he- 
lium burning starts, possibly degenerately in a helium fiash 
(for M ^ 2.25A/0), after which its luminosity decreases. 
For those stars where helium burning starts degenerately, 
the star reaches a higher luminosity on the red giant branch 
than for higher mass stars. This means that although the 
initial luminosity on the giant branch increases with stellar 
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mass, the value at the tip of the giant branch has a minimum 
for stars that no longer start helium burning degnerately 
(M ^ 2.25M0), increasing up to ~ lO^L© for 1.67M0 stars. 
The star's luminosity increases slowly as it moves along the 
Horizontal Branch (HB), which lasts for ~100 Myr. Once 
Helium in the core is exhausted the star swells to become 
an Asymptotic Giant Branch (AGB) star and its luminosity 
increases yet further. Luminosities of evolved A stars can 
reach up to 2 x IO'^'Lq for a 3.8M0 star in this phase. 

The majority of mass loss occurs on the AGB. There is 
a fair degree of uncertainty in mass loss rates and the exact 
m echanism driving the mass loss in AGB stars (see review 
bv lWillso n (2000) ). The stellar evo lution code uses the mass 
loss of iKudritzki fc Reimerd (|l978l ) on the red giant branch 
and horizontal branch: 

Mk = 2 X 10-^-^:^^ Moyr-i (1) 

where M*, R^ and L, are the stellar mass, radius and lumi- 
nosity in solar values. 

On the AGB the formulation of IVassiliadis fc Wood! 

(|l993l ) is used: 

logMvw = -11.4+0.0125 [P-100 max(Af,-2.5, 0.0)] Moyr 

P = min(3.3, -2.07 - 0.9 log A/. + 1.94 log i?,) (3) 

where P is the Mira pulsation period of the star, in days. 
The expansion velocity for the stellar wind is given by: 



vsw 



'13.5 + 0.056P kms" 



(4) 



The mass l oss rate peaks in a superw ind phase, at the 
tip of the AGB (|Vassihadis fc Wood|[l993 ). with 

(5) 



M = 1.36 X 10"^L. M0yr"\ 



These are empirical mass loss rates, fitted to observations of 
RGB and AGB stars. Thermal pulses, with periods ~ lO^yr 
dominate the evolution on the AGB as the star switches 
between helium and hydrogen shell burning. This may lead 
to discrete superwind phases and the multiple shells seen in 
many planetary nebulae. 

The white dwarf core now evolves at constant lumi- 
nosity to higher effective temperature. This luminosity is 
proportional to the white dwarf mass. The hotter radiation 
from the white dwarf core ionizes the expelled gas which 
is seen as a planetary nebula. Once the white dwarf core 
reaches a maximum effective temperature it starts cooling 
at constant radius. The radiat ive cooling o f the white dwarf 
is modeled by Mestel theory (|Meste]||l952l ). with the white 
dwarf luminosity, in solar units, given by 



635MwdZ"- 
[A{t + 0.l)Y 



(6) 



where Mwd is the mass of the white dwarf (Mq), Z is the 
metallicity, t is the cooling age or time since the WD formed, 
in Myr and A is a parameter that is composition dependent. 
In the current models solar metallicity, Z=0.02 and A=15 
for a CO white dwarf is used. Once the white dwarf has 
cooled significantly, crystallization occurs, and the cooling 
rate enters a different phase, and hence in these models a 
changes from 1.18 to 6.48 at an age of 9Gyr, for all spectral 
types. Although prescriptions for white dwarf crystalliza- 
tion ha ve impr o ved s ignificantly sin ce the groundb r eakin g 
work of lMestell (|l952l ). for example iMetcalfe et al] (|2004h . 



However crystalisation onl y occurs when the w hite dwarf 
cool s to ^ 6 000 - 8000 K JMetcalfe et al.ll2004l ) and there- 
fore [Meste]| (|l952h should provide accurate luminosities for 
white dwarfs hotter than this. It is found later in this work 
that debris discs are only detectable around very young, hot 
white dwarfs (see § 15. 2. 311 and there fore differences between 
the cooling theory of iMestell (|l952r ) and more modern pre- 
scriptions are not significant for the current work. 



3 STEADY STATE MODELS FOR 

POPULATION OF DEBRIS DISCS AROUND 
MAIN SEQUENCE A STARS 

3.1 Evolution of an individual disc 

The evolution of an individual disc, due to collisions, can 
be describe d by a simple stead y state model such as that 
presented in lWvatt et al.l l|2007h . These models assume that 
the size distribution of all particles is the same as for an 
infinite coUisional cascade: 

n{D)dD oc D^'^^^dD. (7) 

'- -' This size distribution holds from the largest planetesi- 
mals of size, Dc, down to the smallest dust grains present 
of size, Dmin and assumes that all bodies of all sizes have 
the same strength per unit mass. In the steady state models 
for main sequence stars Dmin is determined by the largest 
grains that are not blown out of the system by radiation 
pressure, the blow-out size (Dbi)- The blow-out size is an 
important property of the system and depends on the stel- 
lar parameters. Assuming spherically symmetric black body 



;rains, of uniform density p in kgm 
BumTelaLlIigTi): 



Dbi = 0.8 



L. \ /2700\ 
Mj \ p ) 



it is given, in pm by 



(8) 



where L, and Af» are the stellar luminosity and mass in 
solar units. 

For such a cascade the majority of the cross-sectional 
area of the disc is found in the smallest particles, but the 
mass of the disc is determined by the largest objects. In these 
models the evolution of the mass in the disc (Mtot in Af®) 
is determined by the coUisional lifetime (tc) of the largest 
objects, with diameter, D^ in km: 

dMtoi -A'f tot 



dt 



tc{Dc)' 



(9) 



The coUis ional lifetime of par ticles of diameter Dc was 
determined in lWvatt et al.l ((20071) and is given by: 



tc 



5.20 X 10" 



^^pr''^Hf)Q7''D. 



5/6, 



Mt 



tAf^'eS/s 



Myr 



(10) 



where r is the radius of the disc, in AU, and dr is the 
width of the disc, taken to be ^, e is the eccentricity of the 
particles, assumed to be equal to their inclination I, QJ, is 
the dispersal threshold for collisions in Jkg~^ and p is the 
density of particles, taken to be 2700kgm~^. So long as the 
only time dependent variable in Eq.[5]is the mass in the disc, 
it can be shown that the mass in the disc evolves as: 



Mtot{t) 



Aftot(O) 

i^ + rj 



(11) 
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For early times, when t « tc, the mass in the disc remains 
approximately constant at its initial value, only turning over 
and falling ofT as Mtot oc i at times t » tc- At late times 
(t in Myr) the mass in the disc tends to a value Mmax that 
is independent of its initial value, 



Mrr^aAt) = 5.2 X 10" 



,pr 



13/3 (-dr 



5/6, 



)Q*^"'D. 



Mt^'\^/H 



Ma 



(12) 



These models are, however, a simplification. It is only 
considered that collisions between the largest bodies change 
the mass in the disc. Only a single value for the dis- 
persal threshold, Qn , is us ed and cratering collisions (e.g 
iKobavashi fc Tanakal ()20ld l) are ignored. However this sim- 
ple procedure allows us to fit the observations. The calcu- 
lated dust luminosities are within an order of magnitude of 
more detailed prescriptions in which collisions between all 
diameter particles are considered, Q^ is a function of di- 
ameter and the siz e distribution is three-phase, for example 
(|Lohne et ai]l2008l). 

In IWvatt et al.l ([20071) the emission properties of the 
disc were calculated using a black body approximation. For 
reasons that are explained in § 14.4.11 a black body approx- 
imati on cannot be used in the current models. The mod- 
els of IWvatt et al.l (|2007|) were updated to incorporate re- 
alistic emission prop erties of grains using the method of 
IWvatt fc Den^ (|2002l '). The prescription of lLi fc Greenberd 



(19971) was used to calculate optical properties for grains 
with a composition of 1/3 silicates, 2/3 organic refractory 
materials and zero porosity, using Mie theory, Rayleigh- 
Gans theory or geometric optics in the appropriate limits. 
The temperature in the disc now depends on the particle's 
diameter, D, in addition to its distance from the star, r in 
AU: 



T{D,r)^ 



(Qabs)T* 
(Qabs) 



T(D.r) 



1/4 



Tbb 



(13) 



where (Qabs)T, and ((5abs)T(D,r) are the particle's absorp- 
tion efficiency averaged over the stellar spectrum and the 
spectrum of a black body radiating at temperature, T, and 
Tbb is the equilibrium temperature of the particle if it were 
a black body, given by: 



Tbb = 278.3 



L 



1/4 



„l/2 



K, 



where L, is the star's luminosity in solar units. The flux 
from the disc, at a wavelength A, is given by: 



Fd.c(A) = 2.98xlO-^^iiP;L jy, 



(14) 



where Mtot is the mass in the disc in Me, Z)min(/im) and 
Dc (km) are the smallest and largest particles in the disc, p 
is the density of particles, in kgm~^, d is the distance from 
the observer to the star, in pc. 



P(A,r) = 



]^bs{X,D)B4\,T{D,r))a{D)dD, (15) 



B,y{\, T{D, r)) is the black body flux, in JySr"^ and a{D)dD 
is the proportion of the total cross-section of the disc found 
in particles with sizes between D and D + dD. 



Parameter 


Orginal value 


New value 


Mmid 


10.0 


7.8 MjB 


Dc 


60.0 


1.9km 


c 


0.05 


0.05 


Qh 


150 Jkg-i 


150Jkg-i 


7 


0.8 


0.8 


^min 


3AU 


6AU 


Tmax 


120AU 


250AU 



3.2 Population of discs around A stars 

The population of debris discs around main sequence A stars 
(spectral type B8 to AO) is relatively well const rained from 
observations by Spit zer at 24/xm and 70/im (jRieke et al.l 
[2,005; Su ct al. 20061'). iRieke et all (|2005l ) searched a sample 
of 76 individual A stars, with ages between and 800Myr, 
for excess with Spitzer at 24/im. This was extende d by look- 
ing at archival IRAS data to a total of 266 stars. ISu et al.l 
(|200y) search a sample of 160 stars for excess at both 24 and 
70^m. Fig. [2]plots the stars that were found to have infrared 
excess in both surveys. Fig. [2] plots the fraction of stars in 
different age bins that are classified to have small, medium 
or large levels of excess. The level of excess decreases with 
stellar age at both wavelengths, but at a faster rate at 24pm. 
The models in IWvatt et all (|2007l ) were fitted to these 
observations, using a population of 10,000 discs, with a 
distribution of initial masses, radii, spectral type, distance 
and ages. The initial masses formed a log normal distri- 
bution centered on Mmid, assu ming the same width as for 
proto-planetary discs, 1.14dex ({Andrews fc Williamsl 120051 ) 
and the initial radii a power law distribution, with the 
number of discs with radius between r and dr given by, 
N{r)dr oc r ''dr, for discs between Tmin and rmax- It was 
assumed that the stars are randomly distributed, evenly in 
spectral type and age and isotropically in distance. In the 
original model the dust properties were assumed to be black- 
body like and the 24/xm and 70/im statistics were fitted by 
adjusting the parameters, Mmid, -De, e, QJ,, 7, Tmin and 

fmax . 

Here we repeated the fit, using the updated formu- 
lation that incorporates the optical properties of realistic 
grains, to find new va l ues fo r this parameter set. As dis- 
cussed in IWvatt et al.l (|2007l ) a degeneracy in the model 
means that these parameters are not uniquely constrained. 
The fit to the data is good so long as both DJ Q*^ e~^'^ 



( Eg. 15 oflWvatt et al.l floO^)) and Mtot,n^^dDc ''^ (Eq.l6 



-1/2 



of IWvatt et al.l (2007])) remain constant. The former can be 
explained by the example of a massive disc composed of 
stronger bodies that evolves on the same timescale as a less 
massive disc of weaker bodies. Alternatively the later can 
be explained by the fact that a massive disc that contains 
more of its mass in larger bodies (i.e. a larger value for Dc) 
has the same dust mass and therefore is equivalently bright 
as a less massive disc where Dc is smaller. Here we chose 
to keep QJj and e unchanged at 150Jkg^^ and 0.05 respec- 
tively, without any loss of generality and performed a fit to 
Mmid, Dc, -y, Tmin and r^^ax- 

In order to fit the observations with a new population 
calculated using the emission properties of realistic grains, 
we want every disc to evolve in the same way as every disc 
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Figure 2. Upp er and middle panel: The fit to tiic observations of 
ISu et al.l 1120061) at 24 (uppe r) and 70/ini (middle), comparable to 
Fig. 2 of lWvatt et al.l |2003). The plots show the fraction of stars 
with flux ratios in different age bins (OOMyr, 30-190Myr, 190- 
400Myr), at 24/im -^^^ = 1-1.25 (diamond: small excess), 1.25-2 
(square: medium excess), >2 (triangle: large excess) and simi- 



larly at lOfivn 



F, 



1-5 (diamond:small excess), 5-20 (square: 



medium excess), >20 (triangle: large excess). Observed values are 
shown with \/~N error bars, whilst model values arc joined with 
dotted, dash and solid lines, for small, medium and large excess. 
Lower: The distribution of planetesimal belt radii. The model 
population is shown with a dashed line, whilst the sub-sample of 
the model population that is detectable at both 24 and 70/jm is 
shown with a dotted line. The observed distribution of radii are 
shown with triangles and y/N error bars. 



in the old population, calculated using a black body ap- 
proximation. The optical properties of realistic grains mean 
that they are hotter and emit less efficiently at longer wave- 
lengths. Therefore in order for the flux ratio ^ "''•; 



^^Sii^ to 

remain the same, disc radii must increase. In fact a good 
fit is achieved by adjusting the disc radii from the values in 
the original model, r2470, to a new value, r = X247o''2470, 
so long as the other parameters are also adjusted accord- 
ingly; thus we keep 7 at 0.8, whilst Tmin and Tmax increase 
by a factor X2470. Then, in order to keep the flux from each 
disc constant, its fractional luminosity (f), the ratio of the 
luminosity of the disc to the luminosity of the star, should 
remain unaltered. Grains are larger than the emission wave- 
length of starlight on the main sequence such that they ab- 
sorb starlight efficiently and / oc f 72°' 1/2 - Therefore we 



need to adjust Mtot and Dc keeping ^'°'^ constant. The 
maximum fractional luminosity that a disc of a given age can 
have due to its coUisional evolution should also remain con- 
stant and thus using Eq. [12] r'''^Dc is also a constant. To- 
gether these mean that Dc and Mmid are altered by X^yg 
and -'i'2470 respectively. 

The conversion between r and r247o was determined for 
realistic grains with the prescribed size distribution around 
main sequence stars and is shown in Fig. [3] There is a func- 
tional dependence of X2470 on r247o and spectral type that 
can be readily understood. For a given spectral type, X2470 
has a minimum at intermediate radii but increases at small 
and large radii. The latter arises because the cooler tem- 
peratures at larger radii mean that the emission at these 
wavelengths is on the Wien side of the black body spec- 
trum. This means that the small increase in the tempera- 
ture of the grains causes a larger increase in the flux ratio at 
these wavelengths and hence the radius inferred, r247o. This 
is compounded by the fact that the temperature of blow-out 
grains falls off m ore slowly with radius than that of black 
body grains (e.g. ISu et al.l (|2009|)). The increase at small 
radii arises because the emission efficiency of realistic grains 
falls off with wavelength, such that the emitted spectrum ap- 
pears steeper (i.e. hotter) than the true grain temperatures. 
This effect is more important where the spectrum is steeper, 
i.e. where discs are hotter. All of these effects are more pro- 
nounced for later spectral type stars because the blow-out 
size is smaller (Eq. [S)) and therefore there is a larger popu- 
lation of grains whose properties depart from black body. 

X2470 was found self-consistently using Fig. [3] to give 
an average value of 2.1. The above discussion suggests that 
an equally good fit to the statistics could be obtained with 
Afmid = 7.8Me and Dc = 1.9km. Indeed, as is shown in 
Fig. [2] this is found to be the case, with a total X24,70,r of 
16.0 compared to 9.8 of lWvatt etahl ((2003). 



4 MODELS OF DEBRIS DISCS AROUND 
POST-MAIN SEQUENCE A STARS 

The steady state models for the evolution of the disc due to 
collisions described in ij3]are extended to include the effects 
of the evolution of the star, described in ij2] The evolution 
of individual discs and changes to their properties, are dis- 
cussed. Individual discs follow different evolutionary paths 
depending on their properties. The plots presented in this 
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tc = 



Figure 3. The ratio of the radius (r) of a disc calculated using 
realistic grains to the radius inferred from the flux at 24^m and 
70/ini, assuming black body emission (r247o), as a function of 
i'2470- 



section are representative, and look at the evolution of a disc 
around a 2.9Mq (equivalent to AO) star, with solar metal- 
licity {Z = 0.02), initial masses in the disc of 1.0, 10.0 and 
100.0 M^ and initial radii of 10, 50 and 100 AU, although 
the population models discussed in Sec. l5.2l use the evolution 
for the appropriate disc and stellar parameters. 



4.1 Radius evolution 

As the star loses mass, the semi-major axes of orbiting bod- 
ies will increase. The star loses the majority of its mass 
in lO^yr at the tip of the AGB. For the majority of discs 
this mass loss timescale is much longer than the orbital 
timescales of particles and therefore the mass loss can be 
considered to be adiabatic. Only for the discs at much larger 
radii do the orbital timescales start to approach the mass 
loss timescales and the assumption of adiabatic mass loss 
may break down, inducing eccentric orbits, for example a 
body initially in a circular orbit at lOOAU will gain an ec- 
centricity of 0.1 if mass loss rates reach W~'^MQyr~^ . 

However these effects are small and are ignored in the 
current models, where the evol ution of the radius (r) d ue to 
adiabatic mass loss is given bv lVillaver fc Livid (l2007J) : 



r{t) = 



r(0)M40) 
M,{t) 



(16) 



Essentially the radius of the disc switches from an initial to 
a final value during the ~ 10^ yr of extreme mass loss, at 
the tip of the AGB, as shown in the upper panel of Fig. [l] 



4.2 Mass in the disc and the collisional lifetime 

As discussed in ijsl the mass in the disc is dominated by 
the largest particles, and hence the timescale on which the 
mass is depleted is dependent on the collisional lifetime of 
the most massive particles, tc (Eg. I10[l. The evolution of tc 
is shown in the middle panel of Fig. [l] Collisions occur most 
frequently (fc is shorter) in the most massive discs, closest to 
the star. The collisional lifetime increases significantly when 
the star loses mass, and is given by: 



tc(0)Mi„t(0)Mr/'(0) 

M,(f)17/3Mtot(i) 



(17) 



Once the mass in the star changes as a function of time, the 
evolution of the mass in the disc is no longer given simply 
by Eq. [O] instead: 



Aftot Ml,{t)Mt^\t) 
oc 



dMtot _ 
dt ^ tc °^ ri3/3(i) 

Using the expression for r(t), 



Ml 



Therefore 



Mu 



Ml'^'^dt 



Mtot(O) 



l + A/tot(0)A'/M, 



17/3 



where 
A'= [1.4x10" 



9.* _ 



5/3, 



dt 



*5/6 



DcQu'r{Qy-"-'M,{Q) 



\13/3i 



(18) 



(19) 



(20) 



''■(21) 



However the mass in the star is approximately constant 
until the AGB, and up until this point Eq. [9] is valid. In 
fact up to the end of the AGB, the evolution of the mass 
in the disc is similar to that on the main sequence in that 
Mtot remains approximately constant for discs with longer 
collisional lifetimes (larger radii) and tends to Mmax fEa. ll2|) 
for discs with shorter collisional lifetimes (small radii), as 
shown in the bottom panel of Fig. [l] Once the star loses 
mass on the AGB, the collisional lifetime of larger objects 
increases so much that it approaches the Hubble time, even 
for close-in discs (see the middle panel of Fig. [4]). Collisional 
evolution is no longer significant for the total disc mass, 
which remains constant throughout the white dwarf phase. 



4.3 Temperature of the disc 

Particles in the disc are heated by stellar radiation. Their 
temperature is a balance between absorption and emission 
and is a function of particle diameter, as shown in the up- 
per panel of Fig. [S] Large grains emit efficiently and act 
like black bodies, whilst medium sized grains (~ /^m) emit 
inefficiently and are therefore hotter than black body. The 
smallest (< /im) grains emit and absorb inefficiently and 
reach a temperature that is independent of grain size. 

To illustrate the evolution of the temperature of par- 
ticles in the disc we show in the lower panel of Fig. [5] the 
change in the temperature of black body particles (appro- 
priate for large grains) and the smallest grains in the disc 
of size Dmin (as calculated in ii4.4|) as the star evolves. This 
follows the luminosity of the star. Along the main sequence 
the temperature of the disc is relatively constant, but it 
increases up to several hundred Kelvin (depending on disc 
radius) 

as the star's luminosity increases on the RGB and AGB. 
The temperature of the disc falls dramatically as the star be- 
comes a white dwarf, mostly because the stellar luminosity 
drops by several orders of magnitude, but also because the 
discs are now further from the star. For this example mass 
star, the only epoch when the temperature of any parti- 
cle in the disc greatly exceeds the black body temperature 
is for white dwarfs, however for later spectral type (lower 
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Figure 4. The evolution of tlie disc radius (upper), collisional 
lifetime (middle) and total disc mass (lower), for different initial 
disc masses, l.OMg (black), lO.OMig (blue), 100. OM^ (red), and 
radii, lOOAU (dashed line), (bottom plot only 50AU (dotted line), 
lOAU (thick line)), around a 2.9M0 star. 



mass) stars the smaller grains may be hotter than black- 
body through all phases of stellar evolution. 



4.4 Smallest particles in the disc 

The size of the smallest particl es (Dmin) in the co llisional 
cascade was determined in the IWvatt et al.l ((20071) models 
and in 33 by radiation pressure (the blow-out size Eq. (S)). 
This is correct for most discs around main sequence stars, 
however there are several other forces that can remove small 
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Figure 5. Upper: The temperature of particles in the disc as a 
function of particle diameter, for a disc initially at 50AU around 
a 2.9Mq star, around a main sequence star (L, = 190Lq), giant 
star (L* = 5OOL0), a horizontal branch star (L* = I5OL0), an 
AGB star (L* = 1.5x10*Lq) or white dwarf (L* = TxIQ-'-^Lq). 
The solid lines show bound grains that are included in the model, 
whereas the dotted lines are unbound grains that are excluded 
from the model. Lower: The evolution of the temperature of black 
body grains (Dc) and the smallest grains in the disc (D^in), for 
discs initially at lOAU (thick line), 50AU (dotted line), lOOAU 
(dashed line) around a 2.9Mq star. 



particles from the disc, including Poynting-Robertson drag, 
stellar wind pressure, stellar wind drag and sublimation. In 
the models presented here, Dmi?i is determined by whichever 
of these removes the largest diameter objects at a given 
epoch. In the following section all five processes are dis- 
cussed and compared such that Dmin can be determined for 
every disc during its evolution. The outcome is summarised 
in Fig. El 

4-4-i Radiation pressure 

Radiation pressure is a radial force which acts in the op- 
posite direction to the star's gravity. The ratio of radiation 
forces to the gravitational forces is given by: 



/3rad ~ 



F,, 



= 0.43- 



M*D 



(22) 



where (Qpr) is the radiation pressure efficiency for 
grains of a given diameter D, in /^m, averaged over the stel- 
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Figure 6. The limiting diameter below which objects are removed 
by cither radiation pressure (Dhl)-, PR-drag (Dpn), stellar wind 
pressure {Dgyy) or stellar wind drag (Dswpr), calculated using 
Eq. [8] Eq. [2S] Eq. [28] and Eq. |30] in a disc initially at lOOAU 
with a mass of lOMg, around an evolving 2.9Mq star. The lower 
panel shows a zoom-in on the AGB. DswPR overestimates the 
diameter bodies that are removed by stellar wind drag as it does 
not take into account the finite AGB lifetime. 



lar spectrum. The dependence of firad on particle diameter, 
for realistic grains, is shown in the upper panel of Fig. [7] at 
different epochs. The peak in Prad occurs at a size compa- 
rable to the peak wavelength in the stellar spectrum. The 
dotted line shows the approximation to /Brad used in this 
work for which (Qpr) = 1. It only deviates from the more 
realistic calculation at the smallest particle sizes and since 
these small particles are generally removed from the disc, 
apart from during the white dwarf phase, this is considered 
a reasonable approximation. 

If /3rad > 1 the radiation forces are greater than the 
gravitational forces and particles are unbound. However 
for small particles produced in collisions radiation pressure 
causes their orbits to differ from that of the parent body, 
resul ting in a hyperbolic orbit if /3rad > 0.5 (jBurns et al.l 
1 1979 ). Radiation pressure also causes bound grains with 
0.1 < /3rad < 0.5 to have eccentricities greater than their 
parents. This causes a perturbat ion to the size distr i butio n 
from that assumed in Eg. [7] ( e.g. IStrubbe fc Chiand (|2006h : 
iThebault fc Augereaul (|2007f )'). but the resulting disc flux 
can be approximated by assuming the eccentricity to be 
constant for all particles and ignoring the effect on the size 



distribution, which is assumed to be truncated at Du cal- 
culated using Eq. |8] 

The middle panel in Fig. [7] shows how the blow-out size 
(Dbi) changes as the star evolves for stars with different 
initial masses. The blow-out size follows the luminosity evo- 
lution of the star, increasing both on the RGB and AGB, 
up to ~cm in size, before decreasing significantly as the star 
becomes a white dwarf. The stellar luminosity of a white 
dwarf is so low that once the star has cooled sufficiently 
{LwD ^ O.ILq), Prad is ucver above 0.5 (as seen in Fig. (Tjl 
and no grains are removed from the system by radiation 
pressure. The change in stellar luminosity or Prad with time 
causes particles to spiral inwards, changing its semi-major 
axis by: 



a(i) = 



a{0){l-Pradm 



(23) 



{1-Prad{t)) ' 

where a(0) and /?rad(0) are the semi-major axis and ratio of 
the radiation pressure to the stellar gravity at the start of 
the white dwarf phase. No particles ever reach the star by 
this process, the maximum fractional change in semi-major 

axis is (1 - Prad{0)). 

This leaves the question of what, if anything, removes 
the smallest particles from collisional cascades in discs 
around white dwarfs, a prob lem which also exists for M- 
dwarfs (jPlavchan et al.l 120091 ). It is possible that magnetic 
effects or interactions with the interstellar medium remove 
the smallest particles in discs around white dwarfs. In this 
work, however, the fate of the smallest particles is left as 
an open question. The bottom panel of Fig. [7] shows that, 
should they exist, particles smaller than Dmin = 10~^m 
contribute negligibly to the total flux in the wavebands con- 
sidered here. This is because, despite the size distribution of 
Eq. [7] meaning that such grains contain the majority of the 
cross-sectional area in the disc, such small grains also have 
extremely low emission efficiencies at longer wavelengths. 
For example even for the extreme case of a disc at lOOAU, 
around an evolved 2.9Mq , with a white dwarf cooling age of 
IMyr and a size distribution that extends down to 10~^''m, 
the contribution of particles less than 10~*m in size to the 
24pm disc flux is only 20%. Thus for practical purposes we 
set Dmin = 10~*^m. 



4-4-2 Poynting-Robertson Drag 

Radiation forces also oppose the velocity of an orbiting dust 
particle, reducing its angular momentum and causing it 
to spiral inwards, by Poynting-Robertson drag (PR-drag), 
changing its radius by order itself on timescales of 



'pD 



Myr. 



(24) 



i.{Qpr) 

Poynting-Robertson drag is only relevant for particle 
sizes for which the PR-drag timescale is significantly shorter 
than the collisional lifetime, since otherwise the particles are 
destroyed by collisions before they have had the opportu- 
nity to migrate. Assuming that the size distribution extends 
down in size indefinitely according to Eq. [7] and that PR- 
drag lifetime varies according to Eq. 1241 both of which are 
valid in the regime where radiation pressure is negligible, it 
is possible to derive a condition for the diameter, DpR, at 
which the collisional cascade is truncated by PR-drag, by 
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Figure 7. The effect of radiation pressure on tfie disc. Calcula- 
tions with (Qpr) = 1 are shown with dotted lines, whilst the solid 
lines show a more realistic calculation. Upper: The ratio of the ra- 
diational to gravitational forces (Prad) for different diameter par- 
ticles in a disc around a main sequence star (L» = 190Lq), giant 
star (L. = 500Lo), AGB star (L. = 1.5 x lO^Le) or white dwarf 
{Lt = 7 X 10~^Lq). The horizontal black line shows firad = 0.5. 
Particles with jS^ad > 0-5 a-re removed from the system by ra- 
diation pressure. The maximum value of firad is less than 0.5 
around a white dwarf, once LwD falls below ~ O.15L0. Middle: 
The change in the blow-out diameter of realistic grains around 
1.67, 2.9 and 3.8 M© or A9, AO and B8 stars. For both plots 
the dotted lines show calculations with (Qpr) = 1, which do not 
vary significantly from the solid lines, which include a more real- 
istic calculation of (Qpr), apart from for small diameter particles. 
Lower: The fraction of the flux per unit log diameter, defined such 
that the area under the curve is 1, for a disc at lOOAU, around an 
evolved 2.9Mq star as a white dwarf with a cooling age of IMyr. 




Figure 8. The maximum diameter object that can be removed 
by stellar wind drag, throughout the disc's evolution, ignoring 
collisions, as a function of disc radius, for a disc around a 1.67Mq 
(A9), 2.9M0 (AO) or 3.8Mq (B8) star. 



comparing the collisional lifetime of the smallest grains to 
their PR-drag lifetime: 



Dr 



= 8.63 X 10" 



^Liif 



dr\2^7/-A^ 



«5/3 



D^ 



M!^^Mlte^o/a 



^m. 



(25) 



In Fig.[6]Dpi?, is compared to Dm (Eq.[8|. Dpa is always 
smaller than Dm for the disc initially at lOOAU with IOMq 
shown, such that collisions and radiation pressure dominate 
over PR-drag which can the refore be ignored, as was previ- 
ously shown in I Wvatj (|2005l ). A similar analysis for discs of 
different mass and radii around different mass stars shows 
that PR-drag can always be ignored except for close-in discs 
or those low in mass. 

PR-drag in discs around white dwarfs is of particular 
interest as a possible mechanism to remove the smallest 
grains. Of potential importance is the fact that objects larger 



than D > 



r(0)2 



can never reach the star due to PR-drag. 



This is because the luminosity of the white dwarf decreases, 
and thus the rate at which objects spiral in decreases with 
time. The critical size is calculated by solving for the rate 
of change of semi-major axis, a, due to PR-drag , for a zero 
eccentricity particle, given by (jBurns et al.lll979l l: 



/£^\ = -^-^""^ (26) 

^dt' STvpc^Da ^ ' 

where L^d is given by Eq. (6] However since collisions still 
occur on faster timescales than PR-drag can act, even for 
the smallest particles present (10~*m, as is seen in the up- 
per panel of Fig. [Sjl , it is not expected that PR-drag has a 
significant effect on debris discs around white dwarfs. 



4.4-3 Stellar wmd pressure 

Mass loss is an important feature of stellar evolution along 
the giant, horizontal and Eisymptotic giant branches, how- 
ever mass loss rates are low and do not affect the disc sig- 
nificantly except towards the end of the AGB. It should be 
noted that there is a relatively large degree of uncertainty in 
the exact mass loss rates of an evolving star, as discussed in 
!j2] The effect of the stellar wind on particles in the disc is 
very similar to that of stellar radiation in that its pressure 
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Figure 9. The change in position of the ice-Hne and silicate sub- 
Hmation radius due to the change in luminosity of a 2.9Mq star, 
Z=0.02, as it evolves. Objects initially outside the ice-line on the 
main sequence, e.g. at lOOAU, could end up inside of the ice-line 
around a giant or AGB star, such that any particles composed 
purely of water ice would sublimate. 



component causes the smallest particles created in collisions 
to have eccentric or unbound orbits and its drag component 
causes them to spiral inwards. Assuming a stationary wind 
model, with a constant wind velocity vsw, the ratio of the 
pressure force due to the stellar wind to the gravitational 
forces is given by; 



Psw = 5.67 X 10 



\oMtVswQs 



(27) 



KLpD ' 

where (Qsw) is the efficiency for momentum transfer from 
the stellar wind, assumed to be 1, iV/» the mass loss rate, in 
Mq yr"""^ and vsw the wind velocity, in kms^'^. 

Just like with radiation pressure grains smaller than 
Dsw or with /3sw > 0.5 would be removed from the system, 
where 



Dsw = 1-13 X 10 



^iMtVsw 
pM, 



(28) 



However, as shown in Fig. [6] Dsw is smaller than Dbi 
throughout the star's evolution, such that removal of grains 
by stellar wind pressure can be ignored. 

4-4-4 Stellar wind drag 

The stellar wind causes particles of diameter, D, in /im, to 
spiral in towards the star on timescales of: 



tsw = 9.4 X lO'^ 



- Dpr'^ 
QswM* 



Myr. 



(29) 



Similarly to for PR-drag these timescales can be com- 
pared to those for collisions Eq.[lO]to derive a condition for 
the diameter below which particles are removed by stellar 
wind drag: 

^2^14/3 d. 2 ,5/3^ 

DswPR = 1M ^yf^- m. 30 

In Fig. |B] Dsw PR is compared to Dm , Dsw and DpR 
throughout the star's evolution. It shows that stellar wind 
drag is only important for the stronger mass loss rates on 



the horizontal and asymptotic giant branches, seen in Fig. [T] 
This is generally true for all disc and star parameters con- 
sidered in this study. 

The strong mass loss rates, however, only act for a rel- 
atively short timescale, shorter than coUisional timescales, 
such that DswPR actually overestimates the size objects 
that are removed by stellar wind drag. Fig. [S] shows the 
maximum diameter particle that can be removed by stellar 
wind drag for a disc around a 2.9Af0 star, given the finite 
AGB lifetime. For the smallest radii discs, since planetesi- 
mals up to Dc (1.9km) are present in our model, almost all 
the mass in the disc is removed by stellar wind drag during 
the AGB phase. However the maximum diameter particle 
that can be removed for large radii discs is not much larger 
than the blow-out size. 

Although the majority of the disc mass at the end of 
the AGB still lies within the main belt, our treatment of the 
effects of stellar mass loss do not include the fate of smaller 
bodies migrating inwards under stellar wind drag. This is 
the main differe nce between our work and numerical sim- 
ulations such as iDong et al.l (|2010f ). iDong et al.l (|2010f ) also 
include the effect of planets on a planetesimal belt, including 
the trapping of planetesimals into mean motion resonances. 
Although our models do not include the flux from small bod- 
ies spiralling in under stellar wind drag in the calculation of 
the disc luminoisities, the amount of material distributed 
between the inner edge of the belt and the star has been 
monitored. This will be discussed further in H5.2.3I in terms 
of the h ot white dwarf discs observed around some stars, for 
example iFarihi et all l|2009l) . 



4-4-5 Sublimation 

As the star evolves to higher luminosities particles heat up 
and may sublimate. For some ideal assumptions, the result- 
ing rate of change o f diameter, D (in m) , is independent of 
the size (|jurall2008h : 



dD 

~dt 



2ao 
P 



To =2SL 
T{t) 



(31) 



where T the temperature, in K, oq = 1.5 x 10^'^kgm~^s~^ 
and To the composition dependent sublimation tempera- 
ture. For pure water ice To = 5, 530 K and for olivine 
To = 65, 300K, meaning that water ices sublimate at ~110K, 
whereas silicates only sublimate at ~ 1,300K. Here we de- 
flne the ice-line and silicate sublimation radius as the radius 
inside of which black bodies have temperatures hotter than 
this. 

The change in the position of the ice-line and silicate 
sublimation radius as the star's luminosity changes is shown 
in Fig. [9] It can be seen that the silicate sublimation radius 
is always smaller than the discs considered in the current 
models, hence the sublimation of silicates can be ignored for 
the population of discs considered. 

Temperatures hot enough for the sublimation of water 
ices, on the other hand, are found in debris discs around 
main sequence stars, for example sublimation of water ice 
is important for comets in our solar system on orbits that 
approach the sun within the ice line of 6 AU. A disc ini- 
tially outside of the ice- line on the main sequence, may be 
inside of it by either the giant or asymptotic giant branches. 
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The sublimation of objects composed entirely of water ice 
would therefore be expected, resulting in significant mass 
loss from objects of all sizes in the disc. Since sublimation 
loss timescales are proportional to diameter, this means that 
smaller objects are removed most rapidly which could trun- 
cate the collisional cascade size distribution at a size larger 
than the blow-out limit. 

However, the behaviour of more realistic objects of 
mixed composition is more complex. Sublimation may not 
proceed at the rate given by Eq. [Sl] indefinitely as wa- 
ter ice below the surface may be protected from sublima- 
tion by the surround ing layers of other non- volatile mate- 
rial (|jura fc XulbOlOl ). As observed for solar system comets, 
sublimation may also lead to the release of small dust grains 
that were originally embedded in the ice, thus increasing the 
number of small grains. Such a process was invoked in the 
models of Ijural (J2004). Thus, although sublimation may 
truncate the size distribution, and so reduce the number of 
small grains, it may also lead to the production of an ex- 
tra population of small grains. Due to this complexity in 
behaviour it is not clear that sublimation cleanly truncates 
the size distribution and it is therefore assumed not to domi- 
nate over other processes in the current models and its effect 
is discussed further in § 15.2.11 



4-4-6 Summary 

The five processes that could potentially remove the smallest 
particles from the disc have been discussed. Fig.[S]provides a 
summary of which processes are relevant as the star evolves. 
Radiation pressure removes the largest particles from the 
disc throughout most of its evolution, apart from on the 
AGE when stellar mass loss rates are high and relatively 
large objects are removed from the collisional cascade by 
stellar wind drag. PR-drag is only relevant for small radii 
discs on the main sequence and giant branch. Uncertainites 
in the outcome of sublimation mean that the models pre- 
sented in this paper assume that discs are unaffected by this 
process; the implications of this assumption are discussedin 
i; 15.2.11 The maximum of Dm, Dswpr, Dpr and Dsw, as 
shown in Fig. [SI was used to determine the cut-off of the 
collisional cascade, Dmin, in our models. For epochs where 
Dmin decreases with time, a time delay would be expected 
before the small grains are replenished by collisions. How- 
ever in these models we assume that this delay is shorter 
than the timescales considered and that the collisional cas- 
cade is instantaneously replenished. For example at the start 
of the white dwarf phase Dmin decreases rapidly, however 
small grains will be replenished quickly by collisions. The 
collisional lifetime for small bodies is short, even though the 
collisional lifetime for the largest bodies in the disc is long 
§113] 



5 OBSERVATIONS 

The preceding section discussed the various processes that 
affect the disc as the star evolves. Combining these pro- 
cesses, together with a knowledge of the change in stellar 
properties as the star evolves (Fig. [T]), the evolution of a 
debris disc and its observable properties on the post-main 
sequence can be modelled. The evolution of an individual 
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Figure 10. The evolution of the total flux from the disc (top 
panel) and the ratio of the flux from the disc to the flux from the 
star (middle panel) at 70/im, as the star evolves. The thick blue 
line in the upper plot is the sensitivity limit of 110/^Jy, whilst in 
the middle plot it shows the calibration limit of Rn^ = 0.1, for 
Spitzer at TOfiva. The star is a 2.9Mq star, with solar metallicity 
(Z = 0.02), at lOpc and the disc has an initial radius of lOOAU. 



disc is, of course, dependent on its initial conditions, char- 
acterized in the current models by its radius, initial mass, 
distance from earth and the spectral type of the star. In this 
section we first consider the evolution of an individual disc, 
with a given set of parameters, and then proceed to discuss 
the evolution of the population of debris discs observed on 
the main sequence around A stars. 

The two quantities of relevance to observations of the 
disc are its flux (Fdisc) and the ratio of the fiux from the disc 
to the flux from the star (7?^). In order for a disc to be de- 
tected by a given instrument, at a given wavelength, its flux 
(fdisc) must be above the sensitivity limit for that instru- 
ment (Fscns), and the ratio of the flux from the disc to the 
flux from the star (R,^) must be above the calibration limit 
{Rv.iim)- The calibration limit is set by the accuracy that 
the stellar flux is known and the quality of the instrumental 
calibration. Here it is assumed that all far-IR measurements 
have the same calibration limit as Spitzer at 70^m (~0.55), 
whilst mid-IR measurements, such as Spitzer at 24pm have 
a calibration limit of ~0.1. Althoug h the instrument cali- 
brati on for Herschel is quoted as 0.1 (|Poglitsch fc WaelkensI 
l2010f ). once the uncertainity in the stellar flux is included 
the limit will be similar to that for Spitzer at 70/im. 
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5.1 Evolution of a lOOAU disc around a 2.9M0 
star at a distance of lOpc 

Fig. 1 101 shows the evolution of -Fdisc and -Rrofim at 70/im, for 
a disc initially at lOOAU, with a mass of 1, 10 or lOOAfg, 
around a 2.9M0 star at lOpc. The blue lines show the sen- 
sitivity and calibration limits for Spitzer at 70/iin, respec- 
tively. The disc is detectable if both iodise and -Rrofim are 
above these limits. From these plots it can be seen that a 
disc of these initial conditions can be detected on the main 
sequence, early on the giant branch and early in the white 
dwarf phase. The exact values of Fdisc and Rrofim , relative to 
the calibration and sensitivity limits, vary significantly for 
discs of difi^erent radius, initial mass, distance from earth 
or around different spectral type stars, however the form 
of these plots, in terms of when Fdisc and -Rtom?™ increase 
or decrease relative to the evolutionary phase of the star, 
remains relatively unchanged. The discussion begins below 
by considering the variation of Fdisc and Rroiim during the 
evolution of a disc with a given set of initial parameters and 
then goes on to consider the changes to this evolution when 
these initial parameters of the disc are varied in § 15.1.11 
§[5T2]and ^ [5X3] 

Along the main sequence the stellar properties change 
only by a small amoun t and the evolutio n of the disc is 
unchanged from that in lWvatt et al.l (|2007l ). The flux from 
the disc falls off with time as collisional evolution depletes 
the mass in the disc. Observations of discs around nearby 
stars with Spitzer at 70/im are in general calibration limited, 
and the example shown at lOOAU is detectable throughout 
the main sequence. 

On the giant branch, the stellar luminosity increases 
by several orders of magnitude (see Fig. [TJ . The increase in 
stellar luminosity heats the disc, however the increase in disc 
flux is small since all the small grains that would dominate 
the emission are removed by radiation pressure (see middle 
panel of Fig. [7]). There is a substantial decrease in R^o^m 
with time along the giant branch, since the increase in stellar 
flux is large, whilst the increase in Fdisc on the giant branch 
is relatively small. The difficulty in observing discs around 
giant stars is therefore the calibration limit, as can be seen in 
Fig. [in] for the example disc for which R^o^m is only greater 
than i?70nm,iim for the first half of the star's giant branch 
evolution. 

As the star moves onto the horizontal branch its lu- 
minosity decreases from the maximum value on the giant 
branch, but remains higher than on the main sequence, 
whilst the stellar temperature remains low (see Fig. [T]). The 
combination of these means that the stellar flux is high and 
Rroiim is small, less than -Rro^m.iim, for the lOM^ example 
disc. This is true for the majority of discs in our population. 

As helium in the core is exhausted, the star swells to 
become an asymptotic giant branch star. It ejects a signifi- 
cant proportion of its mass in a stellar wind and the smallest 
grains are removed by stellar wind drag {Dmin = Dswpr 
Ea. l30p . The stellar luminosity increases and heats the disc 
such that Fdisc remains high, despite the fact that Dmin is 
relatively large. This means that R^ < R^Aim and discs do 
not have an observable excess from a debris disc. However 
AGB stars may be surrounded by expanding circumstellar 
envelopes of material ejected from the star in a stellar wind 



and emission from these dust shells would be significantly 
brighter than a debris disc in the infrared or submm. 

After mass loss ceases, the white dwarf core evolves 
swiftly to higher effective temperature at constant luminos- 
ity, before the stellar luminosity starts to fall as the star 
cools as a white dwarf. As the stellar luminosity decreases, 
Rv increases and it becomes possible to detect emission from 
this example debris disc. For this short evolutionary epoch 
the star is defined as a post-AGB or pre-white dwarf. For 
the purposes of these models we have defined the post-AGB 
phase as the O.lMyr before the start of the white dwarf 
phase. 

There is a sharp drop in stellar luminosity as the stellar 
envelope is ejected and the stellar core is exposed as a white 
dwarf This means that the ratio of the stellar luminosity 
to the disc fiux increases significantly and that observations 
are no longer calibration limited. However the disc fiux falls 
rapidly below the sensitivity limit as the white dwarf cools 
and it is this limit that determines whether a white dwarf 
debris disc is detectable. As discussed in i|4.4.1l even though 
there is no process to remove small dust created in collisions, 
the fiux from these small grains is small and does not make a 
white dwarf debris disc detectable. As can be seen in Fig. llOl 
Fdisc is only greater than Fscns for very young white dwarfs. 



5.1.1 Dependence on disc parameters 

Changes in Fdisc and Ri, with initial disc mass and radius 
are interlinked. Simplistically Fdisc, and thus R^, is propor- 
tional to disc mass and hence discs that are more massive 
are easier to observe. The collisional evolution of material in 
the disc, however, means that there is a dependence of disc 
mass at later times on disc radius, since the collisional life- 
time is shorter for close-in discs than for those further out 
(Eg. nop . For discs that have reached collisional equilibrium 
their mass, at a given age, is independent of their initial 
mass but increases with disc radius fEa. ll2|) . Discs at large 
radii, on the other hand, will not have reached collisional 
equilibrium and so retain their initial mass. This leads to 
the behaviour of Fdisc with radius shown in the upper panel 
of Fig. 111! For close-in discs the disc mass and thus Fdisc 
increases with radius, despite the decrease in disc temper- 
ature. For large radii discs, on the other hand, there is a 
significant variation in Fdisc with initial disc mass and Fdisc 
decreases with radius or disc temperature. This behaviour 
is of particular importance in determining which radii discs 
are the brightest at a given epoch. As can be seen for a 
given individual disc mass, the brightest discs at some later 
epoch are those with intermediate radii at which the largest 
planestessimals are just reaching collisional equilibrium at 
this age. 



5.1.2 Dependence on stellar parameters 

The simplest scaling relation is the distance to the star. 
Rv is unchanged, whilst Fdisc scales inversely with distance 
squared. At large enough distances observations are always 
sensitivity limited. 

In the current models stars with mass between 1.67 and 
3.8 Mq or spectral type A9-B8 are considered. The differ- 
ence between these models that has the greatest effect on 
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the disc is that in stellar luminosity. More luminous stars 
have brighter discs, although this increase in not as large as 
might be expected because the blow-out size also increases 
with stellar luminosity. Thus, the ratio of the disc flux to 
the stellar flux decreases with stellar luminosity, since the 
increase in stellar flux is larger than the increase in disc 
flux. In terms of Fig. [10] this means that, for higher lumi- 
nosity stars, the upper plot is shifted upwards relative to 
the sensitivity limit, whilst the lower plot shifts downwards. 
As was discussed in § [2] the stellar luminosity increases the 
most on the giant branch for lower mass stars and therefore 
-Fdiac and Rv show the greatest variation on the giant and 
asymptotic giant branch for these stars. 

There is also a dependence in disc mass on main se- 
quence lifetime for coUisionally evolving discs. Later spec- 
tral type stars take much longer to evolve and therefore the 
reduction in their disc mass at a given epoch is larger, how- 
ever the difference in disc flux due to this is small compared 
to the difference due to the change in stellar luminosity. 

In terms of detecting discs, this means that where ob- 
servations are calibration limited, disc are more detectable 
around the least luminous stars. This applies to lower mass 
stars on the main sequence, or early on the giant branch. 
On the otherhand when observation are sensitivity limited, 
i.e. around white dwarfs, discs are more detectable around 
the most luminous, or higher mass stars. 



5.1.3 Dependence on wavelength of observations 

The above discussion has focused on observations with 
Spitzer at 70/im. The form of the upper and middle panel 
of Fig. [11] remain relatively unchanged as observations are 
made in different wavelengths however the exact values of 
Fdisc and R^ relative to the sensitivity and calibration lim- 
its vary significantly. The disc flux peaks at approximately 
the peak emission wavelength for a blackbody of the disc 
temperature. The ratio of the disc to stellar flux also has a 
similar variation with wavelength, however it peaks at longer 
wavelengths, since F, falls off more rapidly with wavelength 
than Fdisc- Variations in Fdisc and Rv are larger for shorter 
wavelengths, where the emission is from the Wien region of 
the blackbody spectrum. 

All of this behaviour, means that there will be an op- 
timum wavelength for detecting discs that depends on disc 
temperature, and whether observations are sensitivity or cal- 
ibration limited. When observations are sensitivity limited 
discs are most detectable for the wavelength at which -Fkisc 
is maximum, given by Wien's displacement law for a disc of 
a given temperature, around lOO/xm for young white dwarfs. 
Alternatively when observations are calibration limited, the 
most discs are detectable for the wavelengths at which Ri, 
is maximum, for example on the giant branch this varies 
between 100 and SOO/im. 



5.2 Population models 

Using our models the evolution of a disc, with a given set of 
initial parameters, can be determined. Here we apply these 
models, to the population of discs on the main sequence 
know n from observations of A stars by Spitzer and the mod- 
els of IWvatt et al.l (|2007|) (see §.|3]). These discs are evolved 



from the main sequence through to the white dwarf phase 
and the population of discs around evolved stars is deter- 
mined. The following discussion focuses on giant stars, hori- 
zontal branch, post-AGB stars and white dwarfs. AGB stars 
are not discussed because debris discs are not detectable dur- 
ing this phase, and in any case observations would be com- 
plicated by the presence of material emitted in the stellar 
wind. 

There are many surveys for debris discs with recent and 
current instruments, such as Spitzer and Herschel, as well 
as up-coming instruments such as Alma. Here we calculate 
the percentage of the evolved population that are detectable 
with various instruments. Table [T] shows these percentages 
for young white dwarfs, horizontal branch, giant and main 
sequence stars. 

It is important to note that these percentages only cor- 
respond to the population of evolved A stars, not the en- 
tire population of giants, horizontal branch stars or white 
dwarfs. The number of discs that these percentages corre- 
spon d to can be cal c ulated from the space density of A stars 
from jPhillips et al.l (|2009l ) of 0.0014 ± 0.0001 pc'^ and the 
average main sequence lifetime for A stars, 950Myr, to give 
a density of 3.5 x 10~^pc~'^ for evolved A stars on the gi- 
ant branch, 2.5 x 10~^pc~^ on the horizontal branch and 
1.47 X 10~''pc~^ for white dwarfs with a cooling age of less 
than lOOMyr. It is not possible to tell from observations of 
giant stars whether they are evolved FGK or A stars and it 
is therefore hard to compare the populations, however the 
majority of white dwarfs currently observed are evolved A 
stars and these space densities make a r easonable c o mpar i- 
son with the 80% complete catalogue of ISion et al.l (|2003 ): 
the number of white dwarfs within lOpc less than 100 Myr 
old is predicted to be 0.6 and less than lOOO Myr within 
202C is predicted to be 50, compared to and 30 (|Sion et al.l 
120091 ). This catalogue contains no white dwarfs with cooling 
ages of less than IMyr with 20pc. 



5.2.1 Giant stars 

Early on the giant branch a small set of the evolved popula- 
tion of debris discs have a detectable excess. The following 
discussion defines which discs are detectable, in terms of the 
parameter space specified by initial disc radius, initial disc 
mass, distance to the star, wavelength for observations and 
mass of the star. In order to assess this the disc flux is plot- 
ted as a function of radius in the upper panel of Fig. 1111 As 
discussed in § 15.1.11 this peaks at intermediate radii because 
collisions have depleted the mass in close-in discs, such that 
Fdisc oc Mmax oc r"^'^ (Eq. I12|l , whereas large disc radii re- 
tain their initial masses, and the disc flux falls off with the 
disc temperature or radius. Only discs with R^ > Rum are 
detectable, or those with fluxes above the solid lines in the 
upper panel of Fig. 1111 excluding both small and large radii 
discs. 

This dependence leads to the form of the middle panel 
of Fig. 111! the solid area of which shows the discs that can 
be detected with various instruments. The upper curve is 
the sensitivity limit, whilst the cut-off at low and high radii 
are from the calibration limit. This plot varies with mass 
and age of the star, as well as mass of the disc and wave- 
length for observations. As can be seen for the example disc 
shown, around an evolved 1.67Mq star, of initially lOMg 
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Table 1 . Detection of discs around evolved stars 



at the start of the giant branch, only discs with radii of 
around lOOAU, within ~ 200pc of the sun, are detectable 
with Spitzer, Herschel or Alma at the wavelengths consid- 
ered. As the luminosity of the star increases along the giant 
branch, the distance out to which discs can be detected in- 
creases, however the range of radii of discs with detectable 
excess decreases. This means that the solid (detectable) area 
of an equivalent to the middle panel of Fig. [11] is largest for 
the least luminous giants. The dependence of disc flux on 
wavelength discussed in § 15.1.31 means that the solid (de- 
tectable) area is largest for Herchel PACS at 160/^m. 

A smaller fraction of the population has detectable ex- 
cess on the giant branch than the main sequence, as can be 
seen in Table [l] The significant increase in stellar luminosity, 
compared to the small increase in disc fiux means that fewer 
discs are detectable over the stellar emission {R^, > Rum). 
Spitzer at 24/xm can detect the largest fraction of the popu- 
lation, because observations are calibration limited and the 
calibration limit in the mid-IR is lower than in the far-IR. 
Herschel SPIRE detects the next highest fraction of the pop- 
ulation due to the wavelength dependence of R^ , peaking in 
the sub- mm, as discussed in ij |5.1.3l A sample of stars within 
lOOpc were considered in Table, [l] However, if observations 
with, for example. Alma were made with the intention of 
detecting such discs, a sample that only extended out to 
smaller distances would maximise the rate of detection. 

Our models suggest that around 10% of evolved A stars 
on the giant branch have detectable excess with Spitzer or 
Herschel. This is, however, subject to the unclear effect of 
sublimation on debris discs. Sublimation could have two pos- 



sible effects. Either it removes all small grains, truncating 
the collisional cascade, and thus decreasing the number of 
discs with detectable excess, or it releases a population of 
small silicate grains, increasing the number of giants with 
detectable excess. This makes future observations of giant 
stars with Herschel, in comparision with our models, very 
interesting, as they have the potential to constrain the ef- 
fects of sublimation on discs. 

Our models, however, compare favourabl y with the 
sample of 44 giants brighter than niv — 4.0 mag (jjuralll99d ) 
(see Table [1}, none of which display excess at 60/im (< 3%) 
with IRAS. Infrared excess is, however, found around 300 
of the 40,000 G and K giants in th e Bright Star Catalogu e 
and Mitchigan Spectral Catalogue (|Zuckerman et al.lll995l ). 
although the origin of this emission is not clear. 12 of thes e 
sources are modelled in further detail in (|Kim et aLlbOOlh . 
who suggest that they are more likely to result from spo- 
radic dust ejection or emission from nearby interstellar cir- 
rus rather than black-body grains in a Kuiper-belt disc. In 
order to compare observations of giant stars with our models 
it would be necessary to distinguish between these scenarios, 
potentially with high resolution imaging. 

Another factor that could significantly change the de- 
tectability of discs around giant stars is the radius distri- 
bution of discs in our population. Our models only included 
discs detected with Spitzer at 24/im and 70/^m and therefore 
there is a bias towards small radii discs. This is particu- 
larly relevant, as large radii discs are detectable, particularly 
at longer wavelengths (see middle panel of Fig. \11\) . This 
could be accounted for by extending our models to include 
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Figure 11. Observations of the disc. Upper; The variation in disc 
flux with radius for a 1.67A/q star on the giant branch (2000Myr) 
(black line), horizontal branch (2200Myr) (blue line) and around 
a white dwarf (cooling age of IGyr) (red line), for discs with an 
initial mass of IOMq. The horizontal lines show the calibration 
limits for Spitzer at70^m on the giant branch (black), horizontal 
branch (blue) and white dwarf phase (red), respectively, whilst 
the dotted horizontal line shows the sensitivity limit. Middle: De- 
tection limits for discs of initially lOMgj , around a star of 1.67-Mq , 
at the base of the giant branch (L* = IO.2L0). Excesses can 
be observed for the discs that fall within the solid filled regions 
for Herschel PACS at 160^m, Herschel Spire at 250^m, Alma at 
1mm and Spitzer at lOfim. The dashed filled regions are for discs 
with radii larger than 200AU, not included in the current models, 
that can be detected. A calibration limit of -R„ um = 0.55 is as- 
sumed for all instruments. Lower: same as middle but for a 3.8Mq 
(equivalent to B8) star that has evolved to become a IMyr old 
white dwarf. Note that the disc radii are 3.93 times larger than 
on the main sequence. 



sub-mm observations of debris discs on the main sequence 
or incorporating this radius bias into our modelhng of the 
main sequence population. Including observations of discs 
around main sequence FGK stars would also make our mod- 
els more directly comparable with a sample of giant stars. 
Given these extensions and a technique to distinguish emis- 
sion from a debris discs from stars undergoing sporadic dust 
ejection or emission from nearby interstellar cirrus, it should 
be possible for future observations with Herschel or Alma to 
determine the effect of sublimation on debris discs. 

5.2.2 Horizontal branch stars 



The majority of the discussion in ii l5.2.1l also applies to hor- 
izontal branch stars. Observations are also calibration lim- 
ited, however significantly fewer discs are detectable around 
horizontal branch stars than giant stars, since the stellar 
flux is on average higher, whilst the disc flux remains ap- 
proximately constant. In order to maximise the percentage 
of discs with detectable excess observations should focus on 
low luminosity horizontal branch stars. In Table [T] the per- 
centage of the population of horizontal branch stars within 
20pc and with luminosities lower than IOOLq were calcu- 
lated. The most discs are detectable with Herschel PACs at 
160//m since this is the wavelength at which R,^ is maximum. 

5.2.3 White dwarfs 

As can be seen in Fig. [I0]the disc flux falls off rapidly as the 
star cools during the white dwarf phase and it is therefore 
very hard to detect debris discs around white dwarfs. Obser- 
vations of debris discs around white dwarfs in our baseline 
model are sensitivity limited and only the most massive discs 
around the closest, youngest white dwarfs are detectable. By 
the same analysis as in § 15.2.11 the bottom panel of Fig. 1111 
shows the distance out to which discs of initially lOM® can 
be detected around an evolved 3.8Mq (equivalent to B8) 
star, with a white dwarf cooling age of IMyr. Thus for discs 
in our baseline model (initial radii less than 200AU) around 
white dwarfs that are younger than 1 Myr, the disc flux is 
so low that it is only those that are within a couple of par- 
sees of the Sun that are detectable with Spitzer, Herschel 
or Alma. Even the increased sensitivity of Spica only means 
that discs out to tens of parsec are detectable. 

Similarly to discs around giant stars, it is the large radii 
discs that retain the highest mass at late times, that are 
therefore the most detectable. As can be seen in the upper 
panel of Fig. 1111 -Fkisc peaks at ~200AU for a disc of initially 
IM^ around the IMyr old white dwarf considered. This ra- 
dius increases with initial disc mass or white dwarf cooling 
age. 

There is, however, a balance between young white 
dwarfs being the most luminous and therefore having the 
brightest discs and the low volume density of young white 
dwarfs such that they are more likely to be found at greater 
distances from the sun. Fig.[T2]shows the maximum distance 
out to which discs around white dwarfs can be detected as 
a function of cooling age, for discs at lOOAU with Spitzer at 
70/im, Herschel SPIRE at 250/im and Alma at 450^m. This 
is compared to the distance within which one white dwarf 
of a give n cooling age i s foun d, according to the space den- 
sities of IPhillips et al.l (|2009l ) . The maximum distance out 



16 A. Bonsor et al. 



to which discs can be detected is never significantly greater 
than the distance within which there is one white dwarf 
and it is therefore unlikely that such a system can be ob- 
served. There is an optimum cooling age for detecting white 
dwarf discs, which varies with wavelength, for Spitzer at 
70/im it is ~lMyr, whilst for Herschel SPIRE at 250^m it 
is ~10Myr and Alma at 450/im ~100Myr. As the disc tem- 
perature drops, the disc flux decreases, more rapidly at the 
shorter wavelengths. This means that for a young popula- 
tion of white dwarfs, the best chances of detecting debris 
discs are at the shorter wavelengths of Spitzer or Herschel, 
whilst for a sample that includes older stars Alma would 
be better. However, overall, the best chances of detecting 
such a system are with the longer wavelengths of Herschel 
or Alma. 

Focusing on Spitzer at 70/im, if for some reason our 
models under-predicted the flux from (or mass in) such discs 
by approximately an order of magnitude a disc would be 
most likely to be detected around a white dwarf of less than 
5Myr old at a distance of around 200pc. The only detec- 
tion of excess around a white dwarf that resembles a m ain 
sequence debris disc is the helix nebula (|Su et al.l 120071 ). a 
young white dwarf with a cooling age significantly less than 
5Myr, surrounded by a planetary nebula at 219pc. This fits 
nicely with our models, especially given that alternative ex- 
planations that increase the disc fiux exist, for exa mple the 
trapping of bodies in resonances (|Dong et al.ll2010l ). 

There are very few young white dwarfs close to the sun, 
therefore assuming that our models are correct, the best 
chances of detecting a white dwarf debris disc are to observe 
nearby white dwarfs with Alma. Table [l] shows the percent- 
age of the population of white dwarfs within lOpc with a 
cooling age of less than lOOOMy r. According to the space 
densities of (jPhillips et al.l |2009|) there are only 6 whit e 
dwarfs in this distance, and even less from lSion et al.l l)2009f ). 
and therefore the chances of one of these white dwarfs hav- 
ing a debris disc within the narrow initial radius and initial 
mass range such that it is detectable is slim. Increasing the 
distance limit of the sample does not improve matters as 
at greater distances the discs flux falls below the sensitiv- 
ity limit. Even with the increased sensitivity of Spica the 
chances of observing such a disc around a white dwarf are 
slim. 

These low probabilities for detecting debris discs around 
white dwarfs fits with the fact that Spitzer observations of 
white dwarfs that have only found one white dwarf with 
infrared excess fitted by a disc with a radius of the same 
order of magnitude of main sequence debris discs. There are 
however ~20 observations of hot, dusty discs around white 
dwarfs that are be st fitted by discs o f radii on the orde r 
the solar rad i us e.g . lFarihi et al.l (|2009h . lReach et al.l l|2005h . 
iFarihi et al.l (|2009l ) estimate that 1-3% of white dwarfs with 
cooling ages less than 0.5 Gyr possess hot IR excess. The 
minimum radius of a disc in our population is ~10AU and 
therefore these observations cannot be explained by the discs 
in our population. Material in discs with such a small radius 
will have a very short lifetime and must, therefore, be re- 
plenished. Within the context of the current model we have 
identified a potential source of material for such discs. Stel- 
lar wind drag was included in the current models in as far 
Eis it truncates the coUisional cascade on the AGB. Material 
that leaves the disc will spiral in towards the sun, most of it 
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Figure 12. The solid line shows the maximum distance out to 
which a disc initially at lOOAU, with a mass of lOMgj, around 
an evolved 3.8Mq white dwarf, of a given age can be detected 
with Spitzer at 70/im, Herschel at 250/im and Alma at 450^m, 
whilst the dotted line shows the distance within which there is 
one white dwarf younger than the given age , calcu lated using the 
space density of A stars from IPhillips et al.l l|2009l ) . 



being accreted onto the star during the AGB, however some 
mass will be left between the inner edge of the belt and the 
star, at the end of the AGB. Fig. [13] shows this mass for all 
the discs in our population. The masses in Fig. [13] are sig- 
nificantly higher than the typical dust masses for these ho t 
discs e.g 3.3 x 10~^°Me of GD166-58 (JFarihi et al.ll2008l l. 
and there are even a significant proportion of the population 
with more mass t han the largest su ch disc, GD362, with a 
mass of 0.017Me (jjura et al.li2009[ l. 



However a mechanism 
is still required to move this material in closer to the star. 
This could potentially be scattering by planets inside of the 
disc or the dynamical effects of mass loss on the disc. The 
effects will be considered in more detail in a future piece of 
work. 



5.2.4 Post-AGB or pre-WD stars 



As discussed earlier in Sec 15.1.21 the stellar flux from AGB 
stars is so high that it is hard to detect emission from a 
debris disc. However as the stellar luminosity starts to drop 
just before the start of the white dwarf phase it is possible to 
detect discs around a small proportion of stars (see Table. [l|. 
The analysis is very similar to that for young white dwarf 
stars, discussed in the previous section. The only difference 
is that it is possible to observe discs around somewhat more 
distant stars, however since the post- AGB phase is shorter 
the density of such stars is lower. The chances of detecting 
such a system are therefore slim. Many post- AGB stars have 
IR excess that can be m odelled as a stable Keplerian disc 
(see Ivan Winckell (|2003l l for a review) . However these discs 
are orders of magnitude brighter than a debris disc at this 
epoch and are mos t likely to be connected with binarity 
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Figure 13. A histogram showing the amount of mass left inside 
of the main belt (between r=0 and r = 4ri,eit) at the end of the 
AGB for the population of discs in our models. 



6 CONCLUSIONS 

This work provides a theoretical framework that considers 
all of the effects of stellar evolution on debris discs, firstly 
considering the evolution of an individual disc, before ex- 
tending this to the known population of debris discs around 
main sequence A stars. It is found that debris discs are 
harder to detect around evolved stars than on the main 
sequence. The fraction of discs with detectable excess de- 
creases significantly on the giant branch, yet further on the 
horizontal branch and discs around white dwarfs are very 
hard to detect, although the limitations during this phase 
are different to earlier in the star's evolution. 

The population of discs on the main sequ ence is con- 
strained by Spitzer observations of A stars (JRieke et al.l 
I2OO5I: Su et al.ll2006i) a nd the steady state collisional models 
of IWvatt et al.l ( 20071 ). In this work these models were up- 
dated to include realist ic grains rather tha n the black body 
approximation used in IWvatt et al.l ((20071). This was done 
relatively simplistically by considering that the difference in 
behaviour between realistic and black body grains can be 
explained entirely by an altered disc radius, characterised 
by the ratio X2470 between the radius calculated using real- 
istic grains (r) and the radius calculated using a black body 



approximation (r247o), shown in Fig . [3] This was used to ad- 
just the fit from lWvatt et al.l (|2007l ) and thus the population 
of discs around main sequence A stars was determined. 

In our models debris discs that are observed on the 
main sequence survive the star's evolution, however their 
properties are altered. They evolve collisionally in exactly 
the same manner as on the main sequence, however the 
longer timescales mean that, except for large radii discs, 
their masses are significantly reduced. Discs heat up as the 
stellar luminosity increases on the giant and asymptotic gi- 
ant branches. The increase in disc flux, however, is relatively 
small since small grains are removed by radiation pressure 
and stellar wind drag (on the AGB only). It is shown that 
Poynting-Robertson drag is irrelevant for all discs, includ- 
ing discs around white dwarfs, the only exception being for 
low mass or close-in discs. Adiabatic stellar mass loss means 
that discs around white dwarfs have radii a factor of 2 or 3 
greater than on the main sequence. 



All of these changes in the properties of the disc can 
be put together to determine which discs can be detected. 
In terms of observations of discs around post-main sequence 
stars, the important quantities are the disc flux Fdisc and 
its ratio to the stellar flux, R^, which must, respectively, be 
above the sensitivity and calibration limits of the instrument 
considered. The variation in these are summarised in Fig. 1101 

A smaller fraction of the population can be detected on 
the giant branch than the main sequence. Observations are 
calibration limited. The large increase in stellar fiux at the 
wavelengths considered, compared to the smaller increase in 
disc fiux, means that R^ decreases significantly and discs 
are hard to observe, i?^ decreases with time on the giant 
branch such that only large radii discs around stars early on 
the giant branch have a detectable excess. One limitation 
of our models is the uncertainty in the effect of sublimation 
on the disc, as discussed in H4.4.5I Future observations of 
giant stars with Herschel or Alma, in comparison with our 
models, could potentially constrain the effects of sublimation 
on debris discs. 

Discs around white dwarfs are very faint and thus hard 
to observe. Their luminosity decreases as the stellar lumi- 
nosity falls off with age and the best chances of observing 
a disc are around very young white dwarfs close to the sun, 
however there are very few such objects and thus the chances 
of observing such a system are small. If for some reason our 
models underpredict the fiux from such discs, then the opti- 
mum age and distance for detecting a white dwarf disc with 
Spitzer at 70fim would be at a distance of ~200pc and an age 
of less than 5Myr. This fits nicely with the one detected disc 
around WD2226-210, the young white dw arf at the centr e 
of the helix nebula, at a distance of 219pc (|Su et al.ll2007l ). 

There are however detections of hot dusty discs around 
~20 white dwarfs with radii less than O.OIAU. These are 
a different population from the discs in our models which 
have a minimum radius of 6AU. Within the context of our 
models a potential source of material to replenish such discs 
is material that has left the disc due to stellar wind drag, but 
not yet been accreted onto the star at the end of the AGB. 
However, a mechanism is still required to move this material 
in towards the star. This, and the dynamical effects of mass 
loss on debris discs will be the subject of future work. 
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